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Abstract 

We obtain a generalization of the relativistic diffusion of Schay and 
Dudley for particles with spin. The diffusion equation is a classical version 
of an equation for the Wigner function of an elementary particle. The 
elementary particle is described by a unitary irreducible representation of 
the Poincare group realized in the Hilbert space of wave functions in the 
momentum space. The arbitrariness of the Wigner rotation appears as 
a gauge freedom of the diffusion equation. The spin is described by an 
SU(2) connection of a fiber bundle over the momentum hyperbolic space 
(the mass-shell). Motion in an electromagnetic field, transport equations 
and equilibrium states are discussed. 

1 Introduction 

Approximations leading to a diffusive motion seem to apply to non-relativistic 
as well as relativistic particles. Nevertheless, the problem of a relativistic exten- 
sion of the diffusion theory encounters some difficulties. There have been various 
approaches to a solution of the problem ( see [1] [2] [3] [4] [5] [6] [7] ; for a review 
and further references see [8] [9]). tn our earlier papers [10]- [11] we followed the 
approach of Schay [5] and Dudley [6] who defined a rclatvistic analog of the 
Brownian motion on the phase space as the diffusion preserving the mass-shell 
in the momentum space. We think that one should approach the relativis- 
tic dynamics starting from quantum field theory which describes multiparticle 
interactions in a relativistic way. In [10] [11] we have discussed the approach 
to equilibrium, the friction terms, the transport equations and an interaction 
with electromagnetic field. In another approach, based on a relativistic Wigner 
function [12] [13], transport equations for relativistic particles have been derived 
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from relativistic quantum field theory. These transport equations should have 
a diffusion equation as a classical non-relativistic limit. 

We are interested in a classical diffusion which would approximate relativistic 
motion of a quantum particle with spin. The Wigner function is the quantum 
analog of the phase space distribution. The (quantum) Wigner function in 
general is not positively definite on the phase space although it can be positively 
definite on the position space or momentum space. We have shown in [11] that 
Schay and Dudley definition of the diffusion results from a quantum Wigner 
function which is positive. We have derived the relativistic diffusion of spinless 
particles [11] from a quantum master equation describing an evolution in the 
proper time (a relativistic master equation in proper time is discussed also in 
[14]) 

In this paper we define a dissipativc Lorcntz invariant dynamics of elemen- 
tary particles described by a positively defined probability distribution. This 
is an extension of the diffusion of spinless particles to particles with spin. We 
look for a second order Lorentz invariant differential operator on the mass-shell 
as a generator of the diffusion. Our starting point is an observation that the 
diffusion of Schay [5] and Dudley[6] as well as the diffusion of spinless mass- 
less particles [11] is generated by the operator M^M^ constructed from the 
spin zero representation of generators M^ v of the Lorentz group. An elemen- 
tary particle is described by a unitary irreducible representation of the Poincare 
group. In the Hilbert space of wave functions defined on the momentum space 
the representation is realized by an irreducible representation of the little group 
of Wigner rotations. There is an arbitrariness in the definition of the Wigner 
rotation resulting from the non-uniquness of the boost ((to, 0,0, 0) — > p). As 
a consequence the generators of the Lorentz rotations depend on an arbitrary 
function g(p) G SU(2). The resulting diffusion equation also depends on the 
choice of the boost transformation. We show that the transformation of the 
diffusion equations can be interpreted as a gauge transformation. A gauge in- 
variant description of the diffusion of particles with spin is possible if spin is 
identified with a curvature of an SU (2) vector bundle over the hyperbolic space 
(the mass-shell) . The gauge invariance is used as a guiding rule for a derivation 
of the relativistic diffusion equation in an electromagnetic field. 

In an equivalent description of spin we consider a phase space which is a 
product of the particle phase space with a sphere. The diffusion of a particle 
with spin can be described as a diffusion on the extended phase space (without 
any discrete variables). Then, an evolution of discrete projections a of spin is 
defined by a transition function P T (x,p, a; x',p', a') which is a matrix element 
of the transition function on the extended phase space . 

We discuss transport equations resulting from an elimination of the proper 
time and equilibrium distributions reached at infinite laboratory time. Appli- 
cations of the spin diffusion to a description of relativistic entanglement, deco- 
herence and magnetization are briefly discussed. 
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2 Wigner wave functions of elementary particles 



Let A e SL(2,C) and A be a homomorphism of SL(2,C) onto 50(3,1). We 
consider the one particle states of a particle described by a wave function in the 
momentum space transforming under an irreducible unitary representation of 
the Poincarc group [15] [16] [17] [18] 

U(A,a)^(p) = ^(A,a;p) = exp^ a )^'(^(p))^(A(A- 1 )p), (1) 

where D J is a (2j + l)-dimensional irreducible unitary representation of SU{2). 
The Wigner rotation matrix A(p) in cq.(l) is defined as 

A(p) =a- 1 (p')Aa(p), 

where p' = A(A)p and the boost A(a(p)) of a particle with mass m is defined 
by A(a(p))(mc, 0, 0, 0) = p ( c denotes the velocity of light). 

The wave function in the position space is defined as usual by the Fourier 
transform 

tp(x) = J dpexp(-ipx)ip(p), (2) 

where px = p^x^ = pqXq — px and po = \/p 2 + m 2 c 2 (Greek indices run 
from tO 3, Latin indices from 1 to 3, boldface letters denote three-dimensional 
vectors). The Fourier transform of U (A, a)ip is a non-local function of x because 
A{p) depends on p. 

The representation (1) is unitary in the Hilbert space of functions square 
integrable with respect to dpp^ 1 . There are two Casimir operators character- 
izing irreducible representations of the Poincare group m 2 = P^P^ and the 
Pauli-Lubanski vector 

w^ti/ 1 = -m 2 j(j + 1), 

where 

= \e^ p M»°PP. 
In the neighborhood of 1 we write 

A = cxp CI, 

where fl = iuj^l^, u^" are real parameters and is a representation of the 
Lorentz algebra by 2 x 2 matrices. In vector notation 

fl = iuj rs l rs + uj^^kj = ial — vk. 

Here dj = ej rs uj rs , lj = ^ej rs l rs and Vj — —uioj- We define the generators 
of the unitary representation of the Lorentz group by 

U(A) = exp(iu/^M^). (3) 
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Let us note that there is some arbitrariness in the definition of the generators 
which is a consequence of the non-uniqueness of the boost. We can choose 
a(p) = a(p)g{p) where g(p) is an arbitrary unitary matrix. Then, the Wigner 
rotation is 

A{p) = g-\p')A{p)g{p). (4) 

In the calculations of the generators we differentiate eq.(l) over w^. Then, 
eq.(4) leads to a gauge transformation of the generators. 

First, we fix the gauge choosing (this is the conventional choice [16] [17]) a 
as a Hermitian matrix. In such a case we obtain from eq.(l) 

M^ v =L^ + -£^(p), (5) 

where 



d d 



£ M „ 6 su(2) results from the differentiation of D in eq.(l). If we write Ji = 
\ei mn M mn and Ki = Mi then after differentiation of eq.(l) we obtain 

d_ 

' dp 



Jr = -itrlkPl-^-j: + S r , (8) 



d d 
K r = iPo-g-^ + £rlnPlS n (po + mcy 1 = ipOg-^ + S r, (9) 

where 

[S k , St] = ie Un S n (10) 

is an irreducible representation of su(2). 

Then, the change of the boost (4) leads to a change of the generators (where 
V(g) is the representation of g appearing in eq.(4)) 



V-^OrV -V^LarV (11) 



whereas Mjk does not change because for A e SU (2) the Wigner rotation A{p) 
coincides with A. 

The wave function (2) transforming under an irreducible representation of 
the Poincare group satisfies the wave equation 

( c - 2 a t 2 -v 2 + m 2 c 2 )^ = o. (12) 

We consider random superpositions of states and define the density matrix 

p**>{p,p') = (Mp>PM)) (13) 
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and 

Paa'(x,x')= I dpdp'exp(-ipx + ip'x'){ip (T (p)ip (r/ (p')). (14) 



Clearly, for a free motion 

{cr 2 dj - V 2 + m 2 c 2 )p(x, x') . . 

= (c- 2 d 2 - V' 2 + m 2 c 2 )p{x, x') = 0. [ b) 

We define the Wigner function (matrix) W as the Fourier transform of p 

W aa ,{x,p)=JdkdWj5(p-±k-±V)eMi(k-k')x)p !7!7 ,{k,k'). (16) 

In the position space 

f 1 1 

W CT(7 '(x,p)= / dvexp(-ipv)/w(x+ -v,x- -v). 

Note that from p + = p it follows that 

W+ = W (17) 

and if p > then 



W(p)= J dxVF(x,p) =p(p,p) >0. 



If an observable <J> CT(T ' (p) is a function solely of p then the expectation value 
in quantum mechanics can be expressed in the form 

Tr(p<S>) =Tr j dxdpW(p, x)$(p), (18) 

where the trace on the rhs is over the spin indices. The formula (18) remains 
valid for observables <I>(p,x) on the phase space. However, the Weyl ordering 
must be applied for a definition of a function of non-commuting operators x and 
P 

We can also define the density matrix for a calculation of spin expectation 
values 

Paa> = J dxdpW (7<7 '(x,p). 

The spin matrix in relativistic quantum mechanics is often discussed in relation 
to the EPR experiments [19]. 

The normalization of the probability distribution means 



Tr J dxdpW aa ,(x,p) = 1. 
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We treat an interaction of a particle with an environment as a motion through 
a random medium. We assume that the interaction preserves the Lorentz in- 
variance. In [11] we consider the Stiickelberg proper time formulation [20] [21] 
of quantum mechanics with a dissipation (see also [14]). In such an approach 
the free quantum evolution disturbed solely by a random Lorentz invariant per- 
turbation reads 



2n- 2 {d T p + i[{d% - V 2 + m 2 c 2 ), p\) = M^M^p + pM^M"" - 2M^pM^ 
+m- 2 w ll w^p + pm~ 2 w fl w' J - - 2mT 2 w ll pw il = [M^ v , [M^,p]] +m" 2 [w M , [w^,p\]. 

(19) 

Note that 

[M„„, [M* 1 ", p}} = 3 2 p + pJ 2 - 23p3 - K 2 p - pK 2 + 2KpK. 

From eq.(19) it follows that Trp is preserved by the time evolution with arbitrary 
Hermitian operators M^ v . J and K terms enter eq.(19) with different signs. The 
K terms are in the Lindblad form [22] whereas the J terms are not (there is a 
wrong sign in front of J 2 ). For this reason it is not clear whether the dynamics 
is dissipative. In sec. 3 it will be shown that the dynamics is well-defined and 
dissipativc if the momentum p is on the mass shell H 

2 2 2 2 2 2 2 fnr\\ 

P = Po -Pi ~P2 -Pa = m c ■ ( 20 ) 

In the massless case instead of 2j + 1 states of spin j we have only two states 
of helicity A = j and A = —j. Eqs.(8)-(9) should be transformed to the he- 
licity basis before the limit m — > [23] [24] . The formulae for generators of an 
irreducible representation of the Poincare group in the case m = have been 
discussed in [24] [25] [26] . The different form of the generators results from dif- 
ferent realizations of the little group E{2) of massless particles. We consider the 
formula of Bialynicki-Birula [26] and Shirokov [25] 

d A d 

J r = -ie rlkPl — + -\p\—ln(p a p a ) (21) 

and 

d 

K r = *|Pl^F ~ ^P3^rbPb(PaP a y 1 , (22) 

where a = 1, 2 and A = |p|~ 1 pJ = |p| _1 pS is the helicity equal ±j . 

For later discussion we mention also the formula of Ohnuki ([18] and refer- 
ences quoted there) 

d 

J a = -ieaikPi-7^ + Ap Q (|p| + f> 3 ) -1 , (23) 



d 

Js = -izmPi-Q^ + \ (24) 
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'dp' 

and 



K a =i\p\j— + \e 3ab p b (\p\+p 3 ) 1 (25) 



3 Diffusion equation 

The equation for the Wigner function (16) resulting from eq.(19) reads 

d T W = P »d*W + \n 2 {M^M^W + WM^M^ - IM^WM^) 
= p»d*W + k 2 (j 2 W + WJ 2 - 23WJ - K 2 W - WK 2 + 2KWK 

+m~ 2 cr 2 w ll w tl W + m~ 2 c~ 2 w fi w' J -W - 2m~ 2 c~ 2 w fl Ww' J -^ . 



(27) 



The derivatives over position will have an index x, derivatives without an index 
are over momenta. 

The operators J and K are defined in eqs.(8) and (9). Writing them on the 
rhs in eq.(27) means that the differentiation acts as usual on the function W 
but the multiplication by the matrix S is a multiplication from the right. We 
have for the orbital part 

L^W(x,p) = L^(x)W + L„„(p)W. 

Unfortunately, the action of on W is quite complicated because the spin 
part depends on the momenta. It is determined by a non-local integral kernel. 
We are unable to show whether the evolution (27) preserves the positivity of W 
or not. Eq.(27) simplifies if we restrict ourselves to the momentum distribution 

W{p) = J dxVK(x, p) = p(p, p). (28) 

In such a case the terms L llv (x)W are absent because / dxL /lv (x)W = 0. Then, 
eqs.(27) can be rewritten as differential equations (using eqs.(8)-(9)) 

n- 2 d T W = K~ 2 p»d*W + \m 2 c 2 /\™W - imc( Po + mc)" 1 ^,^,^^, W] 
-Ip 2 (p + mc)- 2 [S, [S,W]] + \{p + mc)- 2 [pS, [pS,W]] = QW, 

(29) 

where 

A™ = Of + dl + &l + {mc)- 2 p J p k &d k + 3(mc)- 2 Pk d k . (30) 

Eq.(29) is equivalent to eqs.(19) and (27) as an equation for W (28)(then the 
term p^d^W = 0). When S = then eq.(29) coincides with the diffusion 
equations of Schay [5] and Dudley [6] [10] [11]. We suggest that eq.(29) defines 
a proper formulation of the relativistic diffusion of particles with spin. Eq.(29) 
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preserves positivity (as will be shown in the subsequent sections) necessary for 
a probabilistic interpretation. In order to preserve the positivity we rejected in 
eq.(29) the L M „(x)W^ terms from the dissipation part in eq.(27) (such a term 
preserves positivity only if x is restricted to x 2 > 0) but we preserved the term 
p^d*W as describing the spatial Hamiltonian evolution . 

W is a Hermitian matrix. We could decompose it into a basis of Hermitian 
matrices. In the case of j = \ 

W = W iQ) + W {k) <r k , (31) 

where a k are the Pauli matrices and W* are real functions. Eq.(29) can be 
rewritten as a set of equations 

K- 2 d T W ( V = K,- 2 p^d^W^ + im 2 c 2 A™H/(°), (32) 
K - 2 d T W = K~ V^W + ±m 2 c 2 A™W 

-mc(p + mc)- 1 (p(d k WW) -p fc VW< fc )) (33) 
-ip 2 (p + mc)- 2 W - \{p + mc)- 2 (pW)p, 

where W = (W^\ W^ 2 \ W^) is a spin current [27] defined in general by 
In the massless case 

Pw - (?(A)V(A'))- (34) 

Here, p and W are 2x2 matrices. The diffusion equation reads (now w^w^ = 
0, we subtract the helicity in eq.(27) instead of the spin;we have calculated 
M^Af" using either eqs.(21)-(22) or (23)-(25),both give eq.(35)) 

K- 2 (d T W -p»d*W) = \A H W (35) 

with 

A H =p jPk d j d k + 3p k d k , (36) 

where k — 1,2,3 and &* — We note that the generator (36) [11] is the 

limit m — > of to 2 c 2 A^ of rcf.[10]. The diffusion equation (35) is the limit 
m — > of the diffusion equation (29) if we set S — > |p| _1 pS. Eq.(35) shows that 
the dissipation does not depend on the polarization. Its inclusion in eqs.(34)- 
(35) is superfluous; the Lorentz invariant dissipation does not depend on the 
polarization in the limit m — ► 0. 
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4 Gauge invariance 

Let us compare the operator on the rhs of eq.(29) with the Laplace-Beltrami 
operator on an SU{2) vector bundle with a connection A over the hyperbolic 
manifold (20) 

A A = g-*{dj - > A ^ur, ,'':,), - iA t ), (37) 

where [10] 

g jl =6 jl +m- 2 c- 2 p>p l , (38) 
<W=<^-f>oW (39) 

and g = m 2 pQ 2 . 

We write the diffusion equation as 

2m- 2 K- 2 c- 2 (d T W - p^dpV) = A A W (40) 

with 

1 -p x S ad (41) 



mc(mc + pa) 
and 

A k g kl Ai = p 2 m- 2 c- 2 ( Pa + mc)- 2 (S ad S ad - p- 2 pS ad pS a<i ), (42) 

where 

S ad $ = [S,$]. 

A similar connection in the momentum space comes from Foldy-Wouthuysen 
transformation in [28] where it is applied to a description of the spin Hall effect. 

The diffusion equation (33) for the spin current W takes the gauge invariant 
form (40) 

2m- 2 K - 2 c- 2 (d T W - p^d^W) = A A W. (43) 

Here, the covariant derivative is dj + Aj with an 5*0(3) (antisymmetric) con- 
nection 

Af = , * : (PbSja - Pa 5 jb ), (44) 

where (AjW) a = Af is a multiplication of the vector W by a matrix A. 
Under the gauge transformation 

Aj = O^AjO + O^djO 

and W — > _1 W. So that eq.(43) stays invariant. 

The diffusion equation (29) defined by M^M^ depends on the choice of the 
generators in the Wigner representation (arbitrariness of the Wigner rotation 
A(p) resulting from the non-uniqueness of the boost). The arbitrariness of the 
generators is expressed by the gauge transformation of the connection A 

Aj = V~ x AjV + iV~ x djV. (45) 
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The gauge (41) is distinguished by the transversality condition 

g^djAk = 0. 

Note that the gauge transformation (45) is different from the one of eq.(ll) 
by a factor of p n . We would get the gauge transformation (11) if we assumed 
the diffusion generator in the form M jM j. The diffusion would have the 
gauge invariance (11) but no relativistic invariance. The subtraction of Mj k Mj k 
changes the gauge transformation. The change comes from the replacement of 
the Podjdj term by m 2 c 2 djdj resulting from the subtraction of Mj k Mj k from 
M jM j. Clearly, this is also necessary for the Lorentz invariance. 

If in eq.(40) we perform the transformation (45) of the connection then 
after the transformation the solution of cq.(40) is rotated into V~ x WtV . A 
calculation of expectation values depends on this gauge rotation. Only gauge 
independent expectation values have a physical meaning. The expectation value 
is expressed by the trace (18). Hence, as local observables $ we should consider 
cither a multiple of an identity matrix or a local function of 

R jk = djA k - d k A 3 + i[Aj,A k ] (46) 

and its covariant derivatives. As non-local observables we could admit polyno- 
mials of the Wilson loop variables 

$c =Y[c^{iA ] dp 3 ), (47) 

c 

where the rhs of eq.(47) denotes a product integral around a closed curve C. 

5 Feynman-Kac-Ito formula 

We derive a probabilistic solution of eq.(29). Let us define the path ordered 
phase factor (a matrix) as a solution of the equation 

dTf = t At o dp l T T? d , (48) 

where p r is the stochastic process on the mass-shell (20) discussed in [10] (the 
circle denotes the Stratonovitch integral [29]). We could describe the stochastic 
process as the diffusion process solving the diffusion equation (32) 

^(x.p)^^ '^,^], 

where the expectation value E[...] is over the diffusion process (x T ,p T ) as has 
been discussed in detail in [10]. The solution T T of eq.(48) is a product integral 
(a unitary matrix) [30] 

S — T 

T? d = ]\exp(iAfodpj,). (49) 

s=0 
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Then, the solution of the matrix equation (29) reads 

W T (x,p) = E{T? d W(x T ,p(T))} - E[T T W{x T ,p T )T-\ (50) 

where in the last formula the matrices are in the irreducible representation (10) 
of su{2), i.e., 

A= 1 pxS. 

mc{mc + po) 

The vector equation (43) has the solution 

W T (x,p) = E[T T W(x T ,p T )}, (51) 

where T T G 0(3) is the product integral with the connection (44). Eqs.(50)-(51) 
give the solutions of the diffusion equations (32)-(33). 

As a consequence we obtain the diamagnetic inequality 

\W r (x,p)\ <E[\W(x T ,p(r))\], 

i.e., the probability distribution of particles with a spin is bounded by the prob- 
ability distribution of spinless particles discussed in [10]. From eq.(50) it follows 
that if the initial value W is a Hermitian positive definite matrix then W T is 
also a Hermitian positive definite matrix, i.e., 



/ 



dpf a (p)W^U>(p) >0. 



Hence, W can be given a probabilistic interpretation which is preserved in time. 

As discussed in sec. 4 the expectation values of S are not gauge invariant. 
For this reason S in not an appropriate variable for a spin. We should use 

Sj = 2m 2 c 2 e : j k iR k i. (52) 
Then, the expectation value of S is 

(Sj) T = 2m 2 <?e Jkl TrfdpE[Tf 1 R kl (p)T T W(p T (p))]- 

It follows that the expectation value is gauge invariant. It can also be seen from 
eq.(50) that (in the absence of an electromagnetic field) the spin evolution is 
just a spin rotation determined by the connection A (41). By differentiation of 
T T in (Sj) T using (48) and the Ito formula [29], we obtain (see also eq.(29)) the 
double commutator 

a r (4) T ~[s,[s,(4) T ]], 

where S is the matrix representation of su(2). Such a dissipative spin dynamics 
is discussed in [31]. 



11 



6 Interaction with an electromagnetic field 



Without external fields the spin is completely described by the Pauli-Lubanski 
vector (in the particle's rest frame w = (0, ms)). It does not change in time. In an 
electromagnetic field the amplitudes do not satisfy the Klein-Gordon equation 
(12). The standard way to include the interaction with an external electromag- 
netic field is to define bispinor amplitudes which satisfy the Dirac equation (the 
space-time derivatives replaced by covariant derivatives) . The bispinor am- 
plitudes describe both a particle and the antiparticle. In an external field the 
division of Dirac wave functions into positive and negative frequencies (electrons 
and positrons) is field-dependent and non-covariant. We need to eliminate the 
negative frequencies. We do it in an expansion in c _1 . Then, the gauge invari- 
ance in momentum space will be our guiding principle for a derivation of the 
final relativistic diffusion equation. We express the lower bispinor by the upper 
one. Then, we expand the field-dependent parts of the Dirac equation in powers 
of c _1 . In the lowest non-trivial order the square of the Dirac operator for a 
particle in an electromagnetic vector potential a ( determining the magnetic 
field B = V x a) and electric field E can be expressed in the form [32] [33] 



D 2 = 8 2 - V k V k + — SB, (53) 
mc 



where 



Vfe = dk - icik + -^^ tkjiEjSi. 

As discussed in [32] [33] we obtain an extended U(l) x SU{2) gauge invariance 
in the position space as a result of the coupling to the spin degrees of freedom. 
The Pauli equation is expressed as the Klein-Gordon equation in a non-Abelian 
gauge field with the symmetry group U(l) x SU{2). 

In the proper time approach to the density matrix evolution we insert the 
Pauli approximation (53) of the square D 2 of the Dirac operator in the commu- 
tator with the density matrix in eq.(19). Then, the Wigner function evolution is 
determined by the motion of position and momenta determined by the equations 

dr m ' 

^ = —F^. (54) 
dr mc 

The equation for the spin comes from the commutator [D 2 ,p] in the master 
equation (19). In the classical non-relativistic limit this equation reads (see a 
discussion in [34] [35] [36] [37] and a generalization to non-Abelian gauge fields in 
[38] [39] [40]) 

dS G 1 

—f- = — £jki(Bk - ekmPrE n )Si. (55) 

clt Tine Attic 
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In the approximation (53)-(55) to eq.(19) the equation (29) for the evolution of 
the Wigner function in an electromagnetic fields takes the form 

n- 2 d T W = K- 2 p»d*W + ±m 2 c 2 A™W - imc( Po + mc)- 1 ^ x V)[S, W] 
-ip 2 (p + mc)- 2 [S, [S,W]} + i(po + mc)- 2 [pS, [pS,W]] 
+^ 2 ^(B - dr c P * E)[S, W] + «- 2 ^(EV + ^B(p x V))W. 

(56) 

We need to add terms of higher orders in c 1 in order to obtain an equation 
which will be gauge invariant in the momentum space and have the correct 
Thomas precession [34] [36] [37] of the spin (it is surprising that the U (1) x 517(2) 
gauge symmetry in the position space contributes to the gauge invariance in the 
momentum space, see also [28]). 

We suggest that the gauge invariant version of eq. (56) with the correct clas- 
sical Thomas precession is 

d T W - p»d*W = ^-m 2 c 2 A A W 

+ 2^(B ~ ^fe)P x E)[S, W] + ^Fj^idj i A 3 )W. 

In eq.(57) the gauge covariant spin S (52) has been introduced. In the lowest 
orders in cr 1 only its non-relativistic version S (as in eq.(56)) appears. The 
term po(po + mc)^ 1 comes from the relativistic theory of spin precession [34]- 
[37] (replacing | of eq.(56)). 



7 Spin diffusion 

In this section we show that eqs.(29)-(33) can be equivalently treated as equa- 
tions for a function W(x,p, n) defined on an extended phase space R 4 xHx 5 2 
(n G 5 2 ). Let g G SU(2) be represented in the form 

g(<j>,6,a)=g 3 (<j>)g2(6)g 3 (a), (58) 

where g 3 corresponds (under the SU(2) — > 0(3) homomorphism) to the rotation 
with respect to the third axis and g 2 corresponds to the rotation around the 
second axis. We can represent a point n G 5 2 as n(</>, 9) = 0(g(4>, 9, a))n 
where n is parallel to the third axis and O G 50(3). Define 

Y ja (0,4>) = Di Q {g(<l>,0,a))- (59) 

For natural j these are the standard spherical functions. Let W be the Wigner 
(2j + 1) x (2j + 1) matrix (16). We define a real function on the extended phase 
space 

W(x,p,n) = Y,Y j A<j>,6)W aa ,(x,p)Y j „,(cj),9). (60) 

era' 
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In the simplest case j = \ it can be checked by elementary calculations that if 
W is of the form (31) then the transformation (60) gives 

W = W {a) +n k W ( - k \ (61) 

where n% = cos#, n\ — sin 6* cos 0, n 2 = sin 9 sin <f). Hence, the expansion (61) in 
rife is equivalent to the expansion (31) in the basis of Pauli matrices a k - In the 
representation (60) the spin operator is represented as a differential operator 

d 

Sk = -ieki r n l -^-p, (62) 

where n 2 = 1. The operators S of eq.(62) are the generators of rotations on the 
unit sphere. They could be represented in the spherical coordinates (6, <j>) in the 
form well-known from the theory of angular momentum in quantum mechanics. 

The gauge transformation (11) is just a rotation of n,i.e., with another real- 
ization of the Wigner rotation, when the solution of the matrix diffusion equation 
(29)transforms as 

W = Vigr'WVig), (63) 

we have 

W'(p,n) = W(p,0(g)n), (64) 

where O(g) is the rotation corresponding to the element g £ SU(2) (it may 
depend on p). 

We can treat eq.(29) as a diffusion equation on the extended phase space. 
The spin evolution (55) can equivalently be described as an evolution of a point 
on the unit sphere 

dn e /T _ 1 , „. 

— = — (B x n (p x E x n . (65 

dr mc 2mc ' 1 y ' 

It leads to the drift 

Y S W = —(b- -J— p x E) (n x V n )VK 

mc V 2mc / 

in the diffusion equation. The spin diffusion part 

A s = p 2 m- 2 c- 2 (p + mc)- 2 (S ad S ad - p- 2 pS ad pS a<i ) (66) 
is expressed as a differential operator 

A s = p 2 m- 2 c- 2 (p + mc)- 2 ((n x V n ) 2 - p~ 2 (p(n x V n )) 2 ) . (67) 
The term 



-imc(p + mc) 1 e r i k pi-^ w [S r , W] 
= mc(p + mc) _1 (n x V n )(p x \7)W 



(68) 
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describes the spin-orbit interaction. Eqs.(29) and (66)- (68) altogether define a 
diffusion on the product of the particle phase space with a unit sphere. 

The complete diffusion equation of the particle with spin in magnetic and 
electric fields equivalent to eq.(56) (with a non-relativistic spin precession) reads 

K- 2 d T W = K- 2 pf*d*W + \m 2 c 2 /\™W + k- 2 -^- c {B - x E)(n x W n )W 

+K~ 2 -^(p x B + mcE)\7W + mc(p + toc)" 1 ^ x V„)(p x V)W 

+ ra g+^ F ((n x V n ) 2 - p- 2 (p(n x V n )) 2 )w. 

(69) 

In sec. 9 we still add a friction drift term to eq.(69)which drives the diffusion to 
an equilibrium. A solution W T (x,p, n) of the diffusion equation (69) starting 
from a non-negative initial condition W remains non-negative. We can relate 
the function W to the matrix W aa i (defined as a Fourier transform of the density 
matrix (16)) using the expansions (31) and (60). On the extended phase space 
the expectation values arc defined as in classical mechanics 

<$(p, n)) T = J dpdnW T {p, n)$(p, n), (70) 

where dn is the invariant measure on the sphere. The gauge invariance means 
the rotation invariance of the integral (70), i.e., it should not depend on 0(p) if 
we make a replacement n — > 0(p)n. 

Eq.(69) is non-trivial even in the non-relativistic limit 

K - 2 d t W = K- 2 pV x W + ±m 2 c 2 AW 

+ fp 2 m- 4 c- 4 ((n x V n ) 2 - p 2 (p(n x V n )) 2 )vK+ ±(n x V n )(p x V))w 

+ ^? B ( n >< Vn)^ - ^fe( P x E)(n x V n )W + ^p x BW + ^EVH^, 

(71) 

where A = d 2 + <9f + df . 

8 Evolution of observables 

The mean values in quantum mechanics are defined in eq.(18) with W as the 
Wigner function. We extend this definition to the diffusion 

($} T =Trj dxdpW T (p, x)®(p, x) = Tr J dxdpW(p, x)® T (p, x). (72) 

If the observable $ does not depend on x then the formula (72) coincides with 
the one of quantum mechanics defining the expectation values by the trace over 
the density matrix (the lhs of eq.(18)). Eq.(72) defines expectation values in the 
theory of diffusion processes. <J> T on the rhs of eq.(72) is the definition of the 
time evolution of $ (the adjoint of the operator W —> W T ). We suggest that 
quantum expectation values are approximated by the expectation values (72) 
taken over the relativistic diffusion process. 
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As discussed in sec. 2 the form of the diffusion equation depends on the 
choice of the Wigner rotation. The expectation values (72) will not depend on 
the choice of the Wigner rotation if the observable <& is covariant under the 
gauge transformation (45). This means that <3? should either be a scalar or 
transform as under the gauge transformation. In the latter case $ can 

be a function of the curvature R (and its covariant derivatives) (46) or the loop 
$c (47). Hence, we must work with variables (a;, p, S) (52) instead of (x, p, S). 
The difference appears only in the relativistic domain because if |p| << mc then 

R jk ~ (2m 2 c 2 )- 1 e jk iSi = (2m 2 c 2 )-%- fe . 

If we choose the boost (4) with the rotation g then the solution of the modified 
diffusion equation will be V(g)~ 1 W T V(g) (where W T is the solution of the dif- 
fusion equation (29) with the Hermitian boost). As an example, according to 

eqs.(45)-(46) the gauge transformed observable $(S) being a local function of 
S is transformed into V _1 <I>(S)V. Hence, the trace in eq.(72) does not depend 
on the choice of the Wigner rotation. 

From eq.(50) it follows that the diffusion with spin is just a unitary rotation 
of the spin followed by the spinless evolution on the phase space, i.e., 

(* T > = j dpdxE[T? d W(p T ,x T )<i>(p,x)} 

= j dxdpE[T-^{p,x)T T W{p T ,x T )\, [<6 > 

where (p T , x T ) is the solution of the same diffusion equations as discussed earlier 
[10] in the model without the spin. 

The evolution (fr T of $ in cq.(72) is determined by the adjoint operator Q* 

K- 2 d T ^ = g*$ = -p»d*$+ im 2 A7$ + lm ^( Po +mc)- 1 e rlk p l [S r ,$] 
-ip 2 (p + mc)- 2 [S, [S, $]] + I(po + mc)- 2 [pS, [pS, $]], 

(74) 

where 

Ag* = df + d\ + flf + {mc)- 2 d j d k p jPk - 3(mc)- 2 d kP k . (75) 

We have two candidates for time in the relativistic diffusion: r (interpreted as 
the proper time ) and x ■ In the non-relativistic limit , when p n ~ mc, xo 
and r enter the solutions of diffusion equations in an additive way. Hence, xo 
is just a shift of r. In the relativistic case in order to obtain a solution of the 
diffusion equation as a function of the laboratory time x we let r — > oo (see 
the discussion in [10]-[11]). In the limit we obtain solutions of eq.(74) which are 
independent of r . The solution of eq.(29)( or more general eq.(57), expressed 
as d T W — QW ) does not depend on r if 

g\V = 0. (76) 

$ does not depend on r if 

5*$ = 0. (77) 
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Eq.(76) is a transport equation in the laboratory time xo which is well-defined 
for xo < ( because it is of the form doW = —AW where —A is a positively 
definite operator). Eq.(77) is a well-defined transport equation for xq > ( 
because it is of the form do<j> = Acfi). If either of eqs.(76)-(77) is satisfied 
then the expectation value (72) does not depend on r. We treat eq.(77) as the 
basic formula determining the evolution of observables measured in experiments. 
Solutions of eqs.(77) and (57) can be related by a random change of time which 
is a generalization of the transformation between the evolution in proper time 
and laboratory time Xq well-known from the relativistic classical mechanics. In 
[11] we have shown that from solutions of eq.(57) (without spin) by a random 
change of time [29] we obtain solutions of the transport eq.(77). In ref. [41] 
the problem is discussed in a mathematically equivalent form but in the inverse 
direction: from the solutions of the transport equation (77) by a random change 
of time the authors derive the solution of the evolution equation in the proper 
time. 



9 An approach to the equilibrium 

We are interested whether the solutions of eqs.(57),(74),(76) or (77) tend to 
a limit (the equilibrium) at r — > oo and xq — > ±oo. In order to achieve the 
equilibrium additional drift terms ( a friction) must be added to the diffusion 
equation. We add a drift 

Y = R~ (78) 

to the diffusion (29). Then 

G Y =G + Y = ^K 2 m 2 c 2 A A + Y. (79) 
Now, the diffusion equation reads 

d T W T = g Y W T . (80) 

The equation for $ is 

a T $ T = g Y *$ T . (si) 

The transport equation (77) reads 

K- 2 p d <£> = K~ 2 pV x $ + ±to 2 c 2 A™*$ 

-imc-^r(p + mc)- 1 e r ikPi[Sr, $] (82) 
-ip 2 (p + mc)- 2 [S, [S, $]] + |(po + mc)- 2 [pS, [pS, $]] - djltf*. 

Eq.(82) could also be written as 

K- 2 Po d a $ = k- 2 P V x $ + im 2 c 2 A(-4)$ - djR?&, (83) 
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where A(— A) means with A — > —A. 

If W is a function (i.e. a multiple of the unit matrix) and we assume that 
$ is a multiple of the unit matrix then the transport and equilibrium equations 
for <3> are the same as in the spinless case [10]. It follows that we obtain the 
Jiittner equilibrium distribution [42] 

<£> EJ = exp(-/3p c) (84) 

if 

Y = {~K 2 0cpoP> + )^j- (85) 
For modified Jiittner distribution [41] [43] 

^ej = Po 1 cxp(-/?p c) (86) 

we have 

Y m = -\k^c PoI P^-. (87) 

The equilibrium distribution (86) does not depend on the spin variables. We 
know magnetized systems which have a non-uniform spin (magnetic moment) 
distribution. We can obtain such a distribution if B = const ^ and E = 0. 
First, let us note that if there is no spin diffusion ( n = ) then 

W B {p, n) = pv 1 cxp(-(3(p Q c + Bn)) (88) 

is the solution of eq.(76) with the generator Q defined by the rhs of eq.(69). The 
spin diffusion [k ^ )also has (88) as the equilibrium distribution if instead of 
the drift (87) we choose 

Y = -i«; 2 /3cpopV + toc/3k 2 (po + mc)^ 1 (n x B)(p x V) 
+/3K 2 ip 2 m~ 2 c~ 2 (p + mc)~ 2 ((n x B)(n x V„) (59) 
-p~ 2 p(n x B)p(n x V„)^. 

The model (88)-(89) may have applications to a description of the time 
evolution of magnetic properties of relativistic systems with spin. 



10 Discussion 

We have discussed an extension of the relativistic diffusion of Schay [5] and 
Dudley [6] to particles with spin. The extension is uniquely determined by the 
rule that similarly as in the spinless case the diffusion should be generated by 
M^ V M^. We obtain a diffusion which is the unique relativistic diffusion in- 
volving both the momenta and the spin. The diffusion defines the Lindblad 
dissipative dynamics of particles with spin. This can be an interesting model 
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for a study of the role of spin in the particle motion in the phase space. The 
relevance of the spin-orbit coupling for dissipation and decoherence has been 
studied before for the non-relativistic velocities [44] [45] . The spin-orbit coupling 
results from the relativistic theory. The diffusive dissipation may be relevant 
to the relativistic description of entaglement and decoherence in the EPR-type 
experiments [19]. The spin diffusion has been studied before in the theory of 
magnetism related to the Landau-Lifshitz equation [46]. However, our relativis- 
tic approach determines the coupling between spin and momenta in the unique 
way. It can be useful in a description of dynamic phenomena in the theory 
of magnetism. The main field of applications should be in the realm of rela- 
tivistic physics. The dissipation resulting from the synchrotron radiation can 
be interpreted as a diffusion in a radiation field . The spin diffusion could be 
compared with experimental results of spin rotation in the synchrotron [47]. 
The relativistic diffusion of particles without spin has already been applied to 
heavy ion collisions [48] . For more precise experiments when the results are not 
averaged over polarizations the description of spin may be relevant. Finally, in 
astrophysics the problem of propagation of neutrinos scattered in the interstel- 
lar medium [49] [50] could be considered as a natural domain of application of 
the diffusion equation. As follows from eq.(35) for the diffusion of light (treated 
as spinlcss in [11]) the spin of the photon does not play any role. 
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